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Abstract. The near-IR He 1 10 830 A transition is a highly sensitive diagnostic for non-LTE effects in astrophysical plasmas. So 
far, non-LTE line-formation computations have failed to quantitatively reproduce observations of this line in the entire range of 
early- A to late-O main sequence stars. It is shown that the non-LTE modelling was insufficient, for the most part, either because 
of inaccurate photoionization cross-sections for the 2s 3 S state or because of neglecting line blocking. New calculations based 
on state-of-the-art atomic data give excellent agreement with observation for the He 1 10 830 A feature, while profiles of the He i 
lines in the visual are retained. 
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1. Introduction 

The 10830A transition (2s 3 S-2p 3 P°) in neutral helium is 
very important for constraining non-LTE effects in the sec- 
ond most abundant element. It involves the metastable state 
of the triplet spin system and is consequently prone to non- 
LTE effects, because the triplets are only loosely coupled via 
collisions to the singlet system. In addition, its location in the 
Rayleigh-Jeans tail of the energy distribution of hot stars gives 
rise to amplified non-LTE effects in this class of objects. 

Since detection of this line in the solar disk spectrum 
(Babcock & Babcock 1934), it has been observed in most kinds 
of stars throughout the Hertzsprung-Russell diagram (HRD). It 
is widely used as an indicator of chromospheric activity in the 
Sun and other cool stars (e.g. Zirin 1982). Towards the hot end 
of the HRD, this feature is observed as a prominent emission 
line, as in O-stars (e.g. Andrillat & Vreux 1979) and in Wolf- 
Rayet stars (e.g. Howarth & Schmutz [T9^2T ). The Hei 10830 A 
transition is also discussed in a galactic and cosmological con- 
text, ranging from local gaseous nebulae to Seyfert galaxies and 
QSOs (e.g. LeVan et al. 1984). The plasma conditions in chro- 
mospheres, stellar winds, and gaseous nebulae lead to rather 
complex physics of the non-LTE line formation (e.g. Avrett et 
al. [T994 Peimbert & Torres-Peimbert IT987b . Therefore it is 
impossible to separate inaccuracies of input atomic data from 
shortcomings in the plasma physics model. Consequently it is 
mandatory to study the transition in a well-understood envi- 
ronment, such as the stellar atmospheres of early-type stars 
with negligible mass-loss. Main sequence stars of spectral 
types early-A to late-O are excellent testbeds. However, non- 
LTE line-formation computations for stars in this region of 



the HRD (Auer & Mihalas [T972TT973l - AM72/73; Dufton & 
McKeith [19801 - DMK80; Leone et al. fl995l- LL P95) have 
at best reproduced observations (Meisel et al. 1982 , Lennon & 
Dufton [13891 - LD89; LLP95) only qualitatively so far. The 
present study resolves the issue. Naturally, the findings will 
have broader implications than for this objective alone. 



2. Model calculations 

Two He i model atoms are used in the present work: i) 
the 'standard' model of Husfeld et al. ( 119891 ). which is a 
slightly improved version of the AM73 model atom, later 
updated with collisional excitation data of Berrington & 
Kingston (1987), and n) an extended 'new' version of the 
Husfeld et al. model. Considering the term structure of the 
new model atom, the difference is that now all states up 
to principal quantum number n — 5 are treated individually 
(previously up to n — 4), with the remainder grouped into 
combined levels for each n in the singlet and triplet spin 
systems. Oscillator strengths are adopted from the NIST 
database (http://physics.nist.gov/cgi-bin/AtData/mainjsd I, 
from Fernley et al. ([1987 FTS87), or from the Coulomb 
approximation, in this order of preference. Photoionization 
cross-sections of FTS87 are employed, with missing data 
assumed to be hydrogenic (Mihalas 1978, p. 99). Effective 
collision strengths for excitation by electron collisions are 
adopted from ab-initio computations of Bray et al. (2000) and 
Sawey & Berrington ( 1993). Additional transitions are treated 
according to Mihalas & Stone (1968), and for the remainder 
of the optically forbidden transitions, the semiempirical Allen 
formula (Allen 119731 is applied. Here collision strengths Q 
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Fig. 1. Comparison of observed equivalent widths for 
He 1 10830 A (Lennon & Dufton IT^l filled circles; Leone 
et al. 119951 open circles) with model predictions: Auer & 
Mihalas JT9731 dotted line), Dufton & McKeith JT98Q1 dashed 
line), and present computations (full lines) for £ = (lower) 
and 8kms _I (upper curve). The abscissa is the reddening-free 
[cj] index. Measurements of LD89 without error bars have un- 
certainties >80mA. Uncertainties for the fast rotators of the 
LLP95 sample can be well above 100 mA. 

varying between 0.001 to 1000 are assumed, based on trends 
derived from examination of the detailed computations of 
Bray et al. and Sawey & Berrington. The lowest D. occur 
for transitions between the ground state and the levels of 
highest excitation energy, while detailed balance is established 
between high-lying energy levels due to the large collision 
rates. Data from Bell et al. (1983 1 have been adopted for 
electron impact ionization of the ground state, along with the 
Seaton ( 1962 > approximation for the other individual levels 
with threshold photoionization cross-sections from FTS87 
where available, or from the hydrogenic approximation. 
Collisional ionization of the packed levels is accounted for 
according to Mihalas & Stone (1968 1. Finally, the theory of 
Barnard et al. dl969> . which is described further in AM72, 
is utilised for a realistic description of line broadening. 
The model for Hen considers all levels/transitions up to 
n = 20. Its details are of no further importance for the present 
investigation. 

The non-LTE line-formation calculations are performed 
in a hybrid approach. Based on hydrostatic, plane-parallel, 
line-blanketed LTE models calculated with the Atlas9 code 
(Kurucz 1993), the non-LTE line formation is solved using 
Detail and Surface (Giddings 1981; Butler & Giddings 1985). 
The coupled radiative transfer and statistical equilibrium equa- 
tions are solved employing an accelerated lambda iteration 
scheme that uses the treatment of Rybicki & Hummer ( 1991 1. 
Line-blocking is realised by considering Kurucz' opacity dis- 
tribution functions. Note that micro turbulence is explicitly ac- 
counted for in the solution of the statistical equilibrium and ra- 
diative transfer by inclusion of an additional term in the width 
of the (depth-dependent) Doppler profiles adopted at that stage 
of the calculations. The grid comprises models for effective 



temperature T e « ranging from 10000 to 35 000 K at a single 
surface gravity value of log g = 4.0 (cgs), which is adequate for 
early-type dwarfs with negligible mass-loss. Solar composition 
is assumed, with the microturbulent velocity £ fixed to two val- 
ues of and 8kms _1 . Theoretical reddening-free Stromgren 
[ci]-indices (where [c\] - c\ - Q.2(b - y)) are determined for 
the Atlas9 model atmospheres using the suite of programs for 
synthetic photometry by Kurucz. The [ci]-index parameterises 
the effective temperature for the comparison with observation. 

The computed He i/n lines at visual wavelengths in several 
of the hotter models (27 500 < T eS < 35 000 K) have been com- 
pared with state-of-the-art non-LTE model atmosphere calcu- 
lations of Lanz & Hubeny (2003 LH03). Excellent agreement 
is found for the He i triplets and He n lines (and also the hydro- 
gen Balmer lines) except for the utmost line cores, indicating 
that non-LTE effects on the atmospheric structure are small at 
line-formation depths and confirming the validity of the hy- 
brid approach. The present computations predict stronger He i 
singlets than LH03, with the discrepancy increasing with de- 
creasing surface gravity. Note that the present modelling gives 
better agreement with the observed visual spectrum of the stan- 
dard star t Sco. The He 1 10 830 A feature is not included in the 
spectrum synthesis of LH03. 

3. Comparison with observation 

New non-LTE line-formation computations have to prove 
their worth by improving the model predictions when con- 
fronted with observation. Spectroscopic investigations of the 
He i 10 830 A transition in early-type main sequence stars 
are scarce, because of the reduced flux of this type of star, 
contamination of the spectra by telluric absorption at these 
wavelengths, and low sensitivity of standard CCDs in the 
J band. Only two sources of measurements of equivalent 
widths Wa from observations with modern instrumentation are 
available from the literature: a sample of 16 slow-rotating 
(v sin i < 30km s" 1 ) B6-09 stars from LD89 and 8 normal Al 
to B2 dwarfs and subgiants from LLP95 (their chemically 
peculiar objects are omitted in the following), almost all of 
them being fast rotators (v sin /> 130km s~'). Both samples 
have y Peg in common, for which satisfactory agreement of the 
equivalent width is found. 

A comparison of equivalent widths from the present non- 
LTE computations with observation is made in Fig. ^ where 
previous model predictions are also summarised, employing 
the same [cj] calibrations as in LD89. As already found 
by LD89, the calculations of AM73 systematically predict 
too strong equivalent widths for Hei 10 830 A, while those 
of DMK80 indicate that this line goes from absorption into 
(strong) emission near spectral type B2. Their calculations 
do not cover effective temperatures below 15 000K. Hybrid 
non-LTE computations based on line-blanketed Atlas9 atmo- 
spheres by LLP95 (not shown) find good agreement in the 
cooler stars, and a qualitatively similar behaviour as DMK80 
predicting slightly less pronounced emission for spectral types 
earlier than B2. Observation, on the other hand, indicates a 
maximum of equivalent widths around spectral type B2, level- 
ling off towards the late-O stars and eventually leading to small 
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Fig. 2. Comparison of departure coefficients bj as a function 
of Rosseland optical depth r ross for the n — 1 and 2 levels of 
He i from computations using the standard (dotted lines) and 
the present (full lines) model atom, and computations for the 
new model atom with the photoionization cross-section for the 
2s 3 S state replaced by that used in the old one (dashed lines). 
Good agreement is found except for the 2s 3 S state, the lower 
level of the He 1 10 830 A transition. Line-formation depths for 
this feature are indicated. The computations are for a stellar 
atmospheric model with T e ff = 30 000 K, log g — 4.0 (cgs), zero 
microturbulence, solar metallicity, and solar helium abundance. 
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Fig. 3. Ratio of the line source function S / to Planck function 
B V (T) at the centre of the 10 830.34 A fine-structure compo- 
nent of the transition, for the same atmospheric parameters as 
in Fig. |2 The same line designations as in Fig.[2]apply. In addi- 
tion, results from a computation using the new model atom, but 
neglecting line blocking (dashed-dotted lines), are also shown. 
In the inset the resulting line profiles are compared after con- 
volution with a Gaussian of FWHM of 20km s _1 . Use of the 
improved photoionization cross-sections turns the absorption 
into emission, while neglect of line blocking effects gives ex- 
cess emission. 



emission, as is the case in t Sco. Consequently, none of the non- 
LTE line-formation computations for Hei 10 830 A from the 
literature reproduce the trend in observed equivalent width of 
dwarf stars over the whole range of spectral types from early-A 
to late-O. 

The new calculations succeed here, finding excellent agree- 
ment with observations, as can be seen in Fig.^ Note that the 
principal error in equivalent width measurements comes from 
continuum placement, such that the uncertainties can be ex- 
pected to be largest in the fast rotators of the LLP95 sample. 
A correct location of the maximum of line strength is indi- 
cated and much weaker emission than in previous cases, reach- 
ing a maximum at T e ff ^31 000 K. The emission is reduced for 
higher microturbulent velocities and levels off towards higher 
temperatures as helium becomes ionized. The observational 
trends imply near-solar helium abundances and low microtur- 
bulent velocities in the cooler stars and enhanced helium and/or 
higher microturbulences for the earlier types (the more mas- 
sive objects). This agrees well with the findings from analy- 
ses of early B-stars in the classical visual wavelength range 
(e.g. Lyubimkov et al. 2004 1 and predictions from recent stellar 
evolution calculations (e.g. Heger & Langer 2000, Meynet & 
Maeder 2000 1. The results are reassuring; however, the reasons 
for the improvements in the predictions need to be understood 
quantitatively. This is investigated in the following. 

4. Discussion 

A comparison of the predictions from several non-LTE calcula- 
tions is made for atmospheric parameters typical of an early B- 



type main sequence star. The computations are based on both- 
the standard and the new model atom, and a version of the latter 
where the photoionization cross-section for the 2s 3 S state is 
replaced by the one used in the old model atom. Departure co- 
efficients bi = rii/n* (the and n* being the non-LTE and LTE 
populations of level i, respectively) for the n = 1 and 2 levels 
of He i are displayed in Fig. [2] Excellent agreement is found, 
except for the 2s 3 S state (and minor differences for the 2p 3 P° 
level). The new model atom predicts level populations by up to 
~10% lower at line-formation depths than the other two mod- 
els, which mutually deviate much less. 

Because of the high sensitivity of the line-source function 
Si to variations of the departure coefficients in the Rayleigh- 
Jeans limit (see e.g. Przybilla & Butler 2004 1, this relatively 
small difference in level populations can result in fundamen- 
tally altered line strengths for Hei 10830 A. While a ratio of 
S i /B v (B v being the Planck function) slightly larger than unity 
in the old model atom (and the modified new one) leads to re- 
duced absorption when compared to detailed balance predic- 
tions, the significantly larger values in the new model give net 
emission (see Fig. [3}- This is primarily because of the consid- 
erable difference in the photoionization cross-sections for the 
metastable state, the FTS87 data being considerably larger, by 
more than 50% at threshold, see Fig. |4] This weakness in the 
AM73 work has already been pointed out by DMK80, but see 
also LD89. Note that the detailed resonance structure of the 
FTS87 cross-sections, which are constrained to a narrow re- 
gion at short wavelengths, are not accounted for in the current 
computations. 
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Fig. 4. Comparison of the photoionization cross-sections for 
the 2s 3 S state of Fernley et al. ( 119871 full line), as used in the 
present He i non-LTE model atom and of Gingerich ( 1964 dot- 
ted line), as used in the AM73 and Husfeld et al. ( 1989 1 models. 
The cross sections differ by more than 50% at threshold. 

The importance of line-blanketed model atmospheres on 
the formation of Hei 10 830 A (despite general systematic ef- 
fects on the atmospheric structure through a steepened temper- 
ature gradient) is less pronounced than indicated by LD89, as 
LLP95 and the present study use the same model atmospheres. 
In fact, the difference between the DMK80 and LLP95 studies, 
on the one hand, and ours on the other, is the allowance for line 
blocking. We therefore computed an additional model neglect- 
ing line blocking as in the older work. Much larger emission 
results from this, enhanced by a factor ~2 in the case discussed 
in Fig. |2] This is facilitated by significantly larger ionization 
of He i by the unblocked UV-radiation, which gives (markedly) 
lower non-LTE populations for the ground state and the n = 2 
levels. In particular, the 2s 3 5 state experiences stronger depop- 
ulation relative to 2p 3 P°, giving rise to the stronger emission in 
He i 10 830 A in the hotter stars. The agreement of the LLP95 
computations with observation for the cooler stars is also ex- 
plained in this picture, as photoionization of He i is inefficient 
there. Note that the profiles of the He i lines in the visual are 
largely insensitive to the effects of line blocking. 

5. Conclusions 

The present study resolves long-standing discrepancies be- 
tween non-LTE computations and observations for the Hei 
10 830 A transition in early-type main sequence stars. Excellent 
quantitative agreement is indicated now on the level of ac- 
curacy of equivalent width measurements from the literature. 
Earlier non-LTE analyses suffered from inaccurate photoion- 
ization cross-sections for the 2s 3 S state and/or neglected line 
blocking. The profiles of the He i lines in the visual are largely 
insensitive to these details. 

Further efforts for constraining non-LTE effects for helium 
are nonetheless required, as the quantitative interpretation of 
a wide range of astronomical objects relies on accurate line- 
formation computations for this element. These should con- 



sider high-S/N and high-resolution observations including the 
complete near-IR range as facilitated by modern telescopes and 
instrumentation. Line-profile information for all available in- 
dicators will allow stringent constraints to be put on the He i 
model atom. 
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